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C/^ ; ABSTRACT 

■ Context. The chemistry in the inner few thousand AU of accreting envelopes around young stellar objects is predicted to vary greatly 

with far-UV and X-ray irradiation by the central star. 
(-H ' Aims. We search for molecular tracers of high-energy irradiation by the protostar in the hot inner envelope. 

Methods. The Submillimeter Array (SMA) has observed the high-mass star forming region AFGL 2591 in lines of CS, SO, HCN, 
I ■ HCN(y2 = l), and HC'^N with 0.6" resolution at 350 GHz probing radial scales of 600-3500 AU for an assumed distance of I kpc. 

O ' The SMA observations are compared with the predictions of a chemical model fitted to previous single-dish observations. 

Results. The CS and SO main peaks are extended in space at the FWHM level, as predicted in the model assuming protostellar X-rays. 
C/3 , However, the main peak sizes are found smaller than modeled by nearly a factor of 2. On the other hand, the lines of CS, HCN, and 

HC'^N, but not SO and HCN(v2 = l), show pedestal emissions at radii < 3500 AU that are not predicted. All lines except SO show 
a secondary peak within the approaching outflow cone. A dip or null in the visibilities caused by a sharp decrease in abundance with 
increasing radius is not observed in CS and only tentatively in SO. 

Conclusions. The emission of protostellar X-rays is supported by the good fit of the modeled SO and CS amplitude visibilities 
including an extended main peak in CS. The broad pedestals can be interpreted by far-UV irradiation in a spherically non-symmetric 
geometry, possibly comprising outflow walls on scales of 3500 - 7000 AU. The extended CS and SO main peaks suggest sulfur 
evaporation near the 100 K temperature radius. The effects of the corresponding abundance jumps may be reduced in visibility plots 
by smoothing due to inhomogeneity at the evaporation radius, varying by +10% or more in different directions. 
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' 1. Introduction phase show a complex spectrum of molecular line emission at 

\.. .... . r- , • , , • (sub)millimeter wavelengths (van Dishoeck 2003). High-Mass 

. ^ X-ray eniission has been reported from high-mass objects of protostellar Objects (HMPOs) are bright in mid-infrared (10-12 

V , the type UCHII and HeAig Ae/Be, representing the late phases ^ ^^^^ ^ ^^^^^.^^ concentration. An open question 

■ of star forniation, with X-ray luminosities in the range 10-"' - ^.^^^^^ important UV photons and X-rays are for regulat- 

■ - ' f '^11 ^ ^^r!".^^ ^^^^^^ ing the chemistry in these earliest phases. In the Hot Core and 

Feigelson 1998; Hofner et al. 2002). Possible origins of the X- j^^pQ j^^^^^ ^ protostellar object is still deeply embedded in 

rays include magnetic activity due to star-disk interaction, ac- ^ collapsing dusty envelope. Thus UV photons and X-rays are 

cretions shocks and shocks in outflows. The high mass makes ^^j^ ^ ^.^^^^^ observable. Their irradi- 

it conceivable that the central object is also an emitter of sub- ^^.^^^ however, may change severely the chemistry of the enve- 

stantial UV radiation. The inner envelope may thus experience (Maloney et al. 1996). 

high-energy irradiation within a few 1000 AU from the central ^ ' , , ... ,t,^ 
object, where processes important in the formation of a star take "ere we focus on AFGL 259 1 , classified as an HMPO or 
place: molecules evaporate, a disk may form, and outflows and ^^^y ^^^^ Core. It is a well studied object and a testbed for high- 
winds are accelerated. With subarcsecond spatial resolution ob- "^^^^ formation. AFGL 2591 is located in a relatively iso- 
servations using (sub)millimeter wavelength interferometers at 1^^^^ I'^gi"" a distance of 0.5 to 2 kpc (as discussed by van der 
high-altitude sites, it has become possible to disentangle these ^ak et al. 1999). A recent estimate puts it at 1 .7 kpc (Schneider et 
phenomena 2006). As this is not confirmed, we will follow the older liter- 
Younger high-mass protostellar objects ai-e particularly in- ^t'^'"^ use 1 kpc in the following for convenience. Numerous 
teresting in terms of their chemisti-y : before the central star is molecular species have been observed in AFGL 259 1 (for a sum- 
formed most molecules are frozen out on dust grains but evap- «^^y ^oty et al. 2002). HCN observed by infrai-ed absorp- 
orate when the young star ignites and starts heating its ambi- ^lon and by submillimeter emission lines reaches a high abun- 
ent envelope. As a result of this, objects in this "Hot Core" '^^"'^^ ^^^ ^ *e '^^"^''^1 '^o'"^' having a radius less than 

175 AU and a temperature of more than 300 K (Lahuis & van 

Send offprint requests to: A.O. Benz, e-mail: benzOastro. Dishoeck 2000; Boonman et al. 2001; Knez et al. 2003). On the 

phys.ethz.ch other hand, CS 5-4 emission has been reported to have a half- 
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power radius of 26 000 AU in single-dish mapping (Carr et al. 
1995). The two extremes indicate the approximate extent of the 
molecular envelope. 

Thermal desorption not only changes molecular abundances 
in the gas phase. It makes these molecules vulnerable to ener- 
getic radiation, converting the predominantly neutral Hot Core 
chemistry into ion-chemistry. Recently, Stauber et al. (2004a, 
2005) have fitted chemical models for AFGL 2591 to column 
densities derived from single-dish line fluxes of 24 molecules 
and included the effects of energetic radiation. The fits improve 
significantly under the assumption of protostellar far-ultraviolet 
(6 - 13.6 eV) emission and substantial X-ray luminosities (6 ■ 
10^' - lO-'^ erg s"'). The spherically symmetric chemical model 
of Stauber et al. (2004a) predicts that far-UV emission reaches 
as far as about 350 AU from the central object. The discovery of 
CO+ and SO+ emission from AFGL 2591 (Stauber et al. 2004b, 
2007) is further evidence for the presence of protostellar far-UV 
emission. X-rays can penetrate farther into the envelope than far- 
UV since the cross section strongly decreases with energy up 
to 20 keV (Morrison & McCammon 1983). Protostellar X-rays 
have noticeable effects in the envelope out to about 5000 AU, 
from where on ionization by cosmic rays dominates. Thus the 
ionization rate changes significantly over the range of a few hun- 
dred AU to a few thousand AU, leaving unique signatures in the 
spatial profile of individual molecular species. 

According to the chemical models of Stauber et al. (2005), 
CS and SO are sensitive to reactions with water plus energetic 
radiation: In the region where T > 100 K, CS and SO are re- 
duced by water in gas phase. However, if this region is irradiated 
by X-rays from the central object, the abundances of CS and SO 
increase greatly. The sulfur-bearing molecules are then enhanced 
by reactions of atomic sulfur - produced by X-ray induced far- 
UV photodissociation - reacting with OH, and the product SO 
reacting with atomic carbon. In a beam with 1" (1000 AU) ra- 
dius, the model abundances of CS and SO are enhanced by more 
than an order of magnitude, if irradiated by a central X-ray lu- 
minosity of some 10^^ erg s The enhancement by protostellar 
far-UV emission is at most a factor of a few for CS and practi- 
cally none for SO. For HCN the result is inverse. Thus, SO is a 
good tracer for X-rays, HCN is a tracer for far-UV. CS is greatly 
enhanced by X-rays, but effectively traces both high-energy ra- 
diations. 

We present new observations of AFGL 2591 with 0.6" (~ 
600 Astronomical Units, AU) spatial resolution. The goal of 
these investigations is to confront the predictions of the model of 
Stauber et al. (2005) reproducing single-dish observations with 
new interferometric measurements that actually probe the rele- 
vant scales. To avoid confusion with low density material plagu- 
ing low-7 data, observations in submillimeter waves are neces- 
sary. Observational details are given in Section 2. The continuum 
observations, briefly summarized in Section 3, are used as a ref- 
erence for the peak location of the molecular lines. The observed 
spatial profiles of molecular species, presented as visibility am- 
plitudes in Section 4, can then be used as a diagnostic for high- 
energy radiations from the YSO. In Section 5, chemical mod- 
els are presented and the relation between geometry and visibil- 
ity amplitudes is studied. Observations and chemical models are 
compared in Section 6 and the results are discussed in Section 
7. Conclusions on the capabilities of chemical modeling under 
high-energy irradiation concerning sulfur chemistry are drawn 
in Section 8. This paper reports first results and will serve as a 
reference for the analysis of future SMA observations and mod- 
eling of AFGL 2591. 



2. Observations 

The SubmilUmeter Array (SMAf] (Ho et al. 2004) on the top 
of Mauna Kea in Hawaii observed AFGL 2591 in the frequency 
bands of 342.6 - 344.6 GHz and 352.6 - 354.6 GHz on May 13, 
2005. The frequency range includes high-7 lines of the sulfur- 
bearing molecules SO and CS, as well as two lines of HCN and 
one of HC'^N (Tab. 1). These lines, sensitive to X-rays or far- 
UV irradiation, allow us to identify protostellar irradiation and 
to separate the effects of X-ray and far-UV emissions. CO^ is 
also in the range, but was not clearly detected. 

The primary frequency resolution was 0.406 MHz (0.5 km 
s ') and uniform across the 2 GHz band. 3C279 was used for 
flux calibration, 3C279 and Jupiter for bandpass calibration. 
Gain calibration was performed using the quasars BL Lac and 
J2015-I-371. The flux calibration uncertainty is of the order of 
10%, but could be occasionally as high as 20%. A total of 5 
elements of the interferometer were used in the extended array 
configuration covering projected baselines from 38 kA to 215 kA 
(32.7 m to 185 m). The FWHM beam size for natural weighting 
is 0.8" X 0.5" with a position angle of 44 degrees. The field 
of view (FWHM beam of single antenna) is 39". The resolu- 
tion in extended array is ideal for the Hot Core envelope from a 
few hundred to a few thousand AU; the interferometric resolu- 
tion allows to probe spatial scales in the emission down to about 
600 AU in radius, and it will miss structures having radii much 
larger than 3500 AU due to missing short baselines. The overall 
telescope time was 5.6 hours in fair weather at PWV x 2 mm. 
The rms noise in the continuum map is 16 mJy/beam. As there 
is no (single dish) total power mode available for the SMA, we 
use previously reported values from the James Clerk Maxwell 
Telescope (JCMT) to estimate the missing flux due to the lack of 
short spacings. 

The data were reduced using the SMA MIR package (Qi 
2005). For image reconstruction the MIRIAD softwar^ (Sault et 
al. 1995) was applied. The continuum emission was sufficiently 
strong to allow for self-calibration. It was subsequently used to 
improve also the image quality of line observations. 



3. Continuum observations 

The observations of the continuum are briefly presented. They 
are used here as the reference for the location of the line emis- 
sion (Fig. 1) and for self-calibration. The results will be dis- 
cussed in a later paper Figure la displays the map of the con- 
tinuum emission after combining both sidebands. A Gaussian 
fit to the emission in the (m, v) plane yields a total flux density 
of 0.662+0.013 Jy, where the error refers to the statistics and 
does not include calibration. The peak flux density in the image 
(Fig. la) is 0.30±0.02 Jy per synthesized beam, identical for the 
uncleaned and cleaned images in natural and uniform weight- 
ing. The peak is located at RA(J2000) = 20''29'"24."880(7) and 
Dec(J2000) = +40°11'19."5(1) as fitted in (u,v) data. The val- 
ues given in parenthesis are the mean error of the last digit re- 
sulting from calibration errors. The position corresponds to peak 
VLA3 by Trinidad et al. (2003) and is also consistent with the 
measurements at 205 GHz by van der Tak et al. (2006). As the 



' The Submillimeter Array is a joint project between the Smithsonian 
Astrophysical Observatory and the Academia Sinica Institute of 
Astronomy and Astrophysics and is funded by the Smithsonian 
Institution and the Academia Sinica. 

' http://sma-www.cfa.harvard.edu/miriadWWW/manuaIs/ 
SMAuguide/smauserhtml/ 
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Fig. 1. Clean maps of AFGL 2591 emissions observed with the SMA. All images display the 348 GHz continuum in grayscale for 
reference. The contours show either the continuum or line emission as indicated in each panel. Contours for the continuum are 3<t, 
6cr, 9cr, etc. For the line data they are 3cr, 5cr, 7cr, etc. Negative isophotes (-3cr) are dashed. The rms noise in the line maps is on the 
average 0.12 Jy/beam. The white asterisk indicates the position of VLA 3 measured at 205 GHz by van der Tak et al. (2006). The 
half-power beam is given in the first image, as well as the length scale assuming a 1 kpc distance. 



peak is not in the phase center, an offset was applied in all sub- 
sequent procedures in (u, v) space. The peak value of the con- 
tinuum emission amounts to a brightness temperature of 10.4 K, 
using Planck's law. The integrated continuum flux at 348 GHz is 
consistent with thermal dust emission and follows the increase 
with frequency reported before in milUmeter waves (van der Tak 
et al. 1999). 



4. Results of line observations 

4.1. Spectra, maps and Gaussian fits 

The spectra of the detected lines at the continuum peak are 
shown in Fig. 2. There are significant dififerences in central ve- 
locities and line widths, but similarities in peak brightness. The 
line centers are consistent with the mean velocity of the enve- 
lope, Vlsr = -5.5 ± 0.2 km s'^ (van der Tak et al. 1999). 
The only exception is HCN, being red shifted by about 2 km 
s~^ in the Gauss fit. A blue-shifted wing extends, however, to 
-13 km s possibly indicating contributions by the approach- 
ing outflow. The widths of the other lines in Fig. 2, particularly 
of HC'^N and SO, are relatively narrow and without compelling 
signatures of outflows or shocks. We note, however, that all lines 
are asymmetric and, except HCN, have their peaks skewed to 
blue velocities. It may be interpreted as the result of a high dust 
column depth in the iimermost part of the envelope or a possible 
disk absorbing the blue-shifted wings. If this is the case, the red 
wing originates from infall. Again the main isotope HCN line is 
an exception and appears to be skewed to the red (Fig. 2). Note 



that CS, SO, and even the main isotope HCN show a dip in the 
spectrum exactly at the systemic velocity. It suggests that these 
line emissions are optically thick. 

Figure 1 also displays maps of the detected lines. Images 
and beam size have been constructed using natural weighting. 
The line data were integrated from -8 to -2 km s~^ except for 
HCN 4-3 where -10 to -0 km s"' was used. These intervals com- 
prise roughly the > 3cr intensities of the hues. The contours 
can be compared with the half-power size of the synthesized 
beam. All main peaks - except possibly HCN(v2=l) - appear 
broader than the beam, consistent with the sizes derived from 
two-dimensional Gaussian fits to the (u, v) data, given in Tab. 1 . 
Thus, most emission is spatially resolved by these interferomet- 
ric observations. Figure 3 shows visibility amplitude plots which 
confirm that most of the flux is resolved out at the observed base- 
lengths. Positional offsets and diameters of the fitted Gaussians 
are also given in Tab. 1 . None of the difli'erences between ma- 
jor and minor axes are significant, and the position angles of the 
source ellipses have errors that do not make the measured values 
statistically significant. Thus it is not possible to determine the 
position angle of the shape of the main peaks. 

All HCN lines and CS show a secondary emission peak 
northwest of the main peak at about the same location. However, 
only in HCN(v2=l) the secondary peak exceeds 5 cr of the back- 
ground rms and is statistically significant. On the average over 
all the HCN and CS Unes, the secondary peak is displaced from 
the main peak by 1180+140 AU at an angle of 21.5° from the 
RA-axis. This angle may be compared to the opening angle of 
the outflow as seen in near-infrared, using the speckle interfer- 
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Table 1. Observed molecules and properties in the area defined by the SMA synthesized beam centered on the maximum position 
of the map: peak brightness in [Jy/beam], width of Gaussian line shape fitted in velocity and shift Vlsr with relation to the local 
standard of rest, positional offset of the main peak from the continuum peak, and source size. The values in parentheses indicate the 
error in units of the last given decimal(s). ^ Peak brightness of synthesized beam at position of peak continuum, using a Gaussian 
fit in velocity. 
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Fig. 2. Spectra of AFGL 2591 line emissions detected with the 
SMA at the continuum peak. The brightness temperature is aver- 
aged in an area equivalent to the synthesized beam. The dashed 
vertical line indicates the systemic velocity of AFGL 259 L 



ometry image of Preibisch et al. (2003). At a distance of 10 - 
20", we measure a half opening angle of 30±2° relative to the 
RA-axis. Thus, in projection, the secondary emission peak lies 
within the approaching outflow cone. 

The continuum and SO line do not show enhanced emission 
at the position of the secondary peak. The main peak of the SO 
emission however is particularly extended and deviates consid- 



erably from a circular shape at the 5 cr level. The apparent devi- 
ations of SO are to the west (approaching outflow), but possibly 
also to the south. These observations suggest that such deviations 
exist, but more observations are needed to confirm and analyze 
them in space and velocity. SO is the only line whose emission 
is marginally displaced in space from the position of the contin- 
uum emission in the southern direction. 

As the line emissions have significant spatial dimension in 
SMA images. Tab. 2 lists the fluxes contained in the synthe- 
sized SMA beam centered at the peak of the continuum emis- 
sion and compares it to the total flux. The central beam refers 
to the brightness of the central region with a radius of approxi- 
mately 300 AU integrated over the line from -15 km s ' to +4 
km s"' . The total Une fluxes observed by the SMA are integrated 
in a square of 3" x 3". We found this the largest feasible area, 
beyond which the noise contribution dominates. The flux of the 
central SMA beam is a smaller fraction of the total SMA flux, 
the larger the source. The ratios are different for the observed 
molecules: low for CS and SO, higher for HCN 4-3, HCN(v2= 1) 
and HC'^N. The small ratios of CS and SO suggest that the main 
CS and SO emissions have an effectively larger size, and bright- 
ness is well extending into the 600 - 3500 AU range. This can 
be noticed also directly in Fig. 1, comparing CS and SO with the 
HCN fines. 

Table 2 furthermore juxtaposes the spectrally integrated 
fluxes within the central SMA beam and the JCMT beam, prob- 
ing different spatial scales. SMA fluxes are between 3 and 40 
times smaller than the JCMT values measured by a beam that is 
530 times larger in area. If the emission came from a peak unre- 
solved in both beams, the fluxes would be the same. The bright- 
ness ratio of HCN(v2=l) is lowest, suggesting that the source 
of this line is not much larger than the SMA beam, as noticed 
before (Fig. 1 and Tab. 1). The JCMT fluxes for the main HCN 
isotope and CS are more than 30 times larger, indicating that 
their main contribution originates from the envelope at several 
thousand AU, the central peak being probably limited by optical 
thickness. 

The widths (FWHM) of the velocity distributions are also 
given in Tab. 1. In the JCMT field of view (14" diameter), tiie 
line widths are relatively uniform in the range of 4.3+1.0 km 
s~^ The values measured by SMA, referring to a region more 
than an order of magnitude closer to the YSO are around 6 km 
s"^ HCN(v2 - I) is noteworthy as it changes the least. As 
HCN(v2 = I) is very compact in the SMA observations, this 
can again be interpreted in terms of the emission being domi- 
nated by a centrally peaked component. It adds to the claim that 
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Table 2. Observed properties of line observations. The line flux in the central SMA synthesized beam of 0.68" x 0.49" and total flux 
in the SMA image are compared with the results from van der Tak et al. (1999, 2003), who used the 14" beam width of the JCMT. 
The corresponding values of the X-ray model (Section 5.1) are also given for 0.5" and 14", respectively. The visibility amplitudes 
in Jy km s~^ in the intervals 30 - 40 kA and 200 - 220 kA are listed as observed and modeled. 



for this line the JCMT is picking up mostly the same flux seen 
by the SMA (Tab. 2). We note for comparison that the free-fall 
velocity from infinity to 7000 AU, relevant for the JCMT beam, 
amounts to 2.6 km s"'. At 500 AU corresponding to the SMA 
synthesized beam, it is 7.5 km s neglecting the mass of the 
iimer envelope and disk. 

4.2. Brightness distributions (visibilities) 

The visibiUty amplitudes vs. (m,v) distance, representing the 
brightness distributions in space are presented in Fig. 3. The 
(m, v) distance, XkA, corresponds to the baseline length b in units 
of the observed wavelength A projected to the plane of the sky 
at the peak position. It is related to the angular structure by the 
Fourier transformation of the image. The angular fringe separa- 
tion A0 is given by 



AO = 



bsinO 



206.26 



[arcsec]. 



(1) 



When discussing spatial structures, it is more convenient to use 
angular radii r, where 2r = AO, and convert to AU (1" = 1000 
AU for our assumed distance to AFGL 2591). Each baseline 
yields an average value in one dimension. In spherical symme- 
try, the average source diameter is J = 0.5nr. Thus the relevant 
radius in the model is related to the one-dimensional resolution 
by 



2A6' 
n 



131.31 

XkA 



[1000 AU]. 



(2) 



The visibilities at the longest and shortest baselines yield in- 
formation on the smallest observable spatial dimension. At the 
longest baselines (> 200 kA, see Fig. 3), all observed lines, ex- 
cept for SO and possibly the main isotope HCN, show visibil- 
ity ampUtudes exceeding the noise. This indicates that the map 
contains a peak, most likely in the center, that is partially un- 
resolved (< 600 AU in radius). Table 2 gives the average value 
in the range 200 - 220 kA, equivalent to the unresolved flux. As 
HCN(v2=l) and HC^^N have an unresolved component, its rel- 
atively low value in the main HCN isotope is likely an optical 
depth effect. 

5. Models 

The visibility amplitudes of lines contain information on the 
radial variations of molecular abundances. In this section, the 



tools and models are presented that will be used in the follow- 
ing section to compare theory and observations. First, the results 
of chemical modeling are shown, assuming a spherically sym- 
metric density and temperature envelope model derived from 
observations of AFGL 2591. The chemical models include ice 
evaporation and assume different levels of far-UV and X-ray ir- 
radiation. The models predict visibiUty ampUtudes that will be 
compared to observations in the following section. 

Models of the fractional abundances of CS, SO and HCN 
are shown in Fig. 4a-c. Molecular abundances are calculated in a 
time- variable chemical network including 395 molecular species 
and 3865 reactions (Doty et al. 2002). The input density model 
(Fig. 4d) is based on a 1/r power-law distribution derived by 
van der Tak et al. (1999), and the input temperature model is 
from Doty et al. (2002). Abundances are determined in time and 
radial distance. A first model, which we call the 'no-X' model 
in the following, fits best 24 column densities derived from ob- 
servations under the assumption of no internal X-ray or far-UV 
radiation. The deviations reach a minimum for a chemical age 
of 5-lOVr (Doty et al. 2002). Stauber et al. (2005) expanded 
the model to allow for irradiation by X-rays and far-UV and 
searched for the best fit for the same molecular observations, 
allowing for irradiation by a central source. The fitted param- 
eters are X-ray luminosity, the temperature of the X-ray emit- 
ting plasma, far-UV irradiation at 200 AU, and age. The result, 
which we call 'X-ray model' although it also includes far-UV 
irradiation, requires Lj = 8 ■ lO"'' erg s"\ 10^ K, Go = 10, and 
5-10'*yr, respectively. Go is the scaling factor of the average inter- 
stellar far-UV radiation flux (Go=l corresponds to 1.6- 10^^ erg 
cm"^s"' at 6< hv <13.6 eV). Comparing modeled and molecu- 
lar column densities observed by single dishes, the X-ray model 
has a significantly smaller deviation than the no-X model 
(Stauber et al. 2005). 

The dominant form of sulfur in the grain mantles depends 
on the gas density of the pre-stellar phase. Low density favors 
hydrogenation and H2S (Charnley 1997), and high density leads 
to oxygenation (OCS, van der Tak et al. 2000). Laboratory re- 
sults by ColUngs et al. (2004) suggest that the sulfur-bearing 
molecules are small enough to diffuse through the structure of 
porous amorphous water ice and desorb at intermediate temper- 
atures between 100 K (H2O) and 20 K (CO) under interstellar 
conditions. In the following we refer to the location of H2S and 
OCS evaporation as the 'sulfur evaporation radius' . Within this 
radius, sulfur resides in gas-phase molecules, rapidly undergo- 
ing reactions that create CS, SO, SO2 and other sulfur-bearing 
molecules. Outside of the sulfur evaporation radius, sulfur is as- 
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Fig. 3. Visibility amplitudes of AFGL 2591 in the observed line emissions vs. projected baseline distance in units of 10^ A. The 

observed values are shown with error bars. They were integrated from -9.5 to -1.5 km s~' except for HCN 4-3 where -10.5 to 
-0.5 km s"' was used. The dotted histogram indicates the expected noise level assuming no signal. The solid curve represents the 
model (introduced in Section 5), assuming protostellar far-UV and X-ray emission. The model without X-rays is shown dashed. 
The single-dish values measured at JCMT are beyond the scale except for HCN(v2=l), where it is indicated by an asterisk at zero 
(m, v) distance. The visibility amplitudes of the continuum emission (integrated in space, Fig. 3f) are given in flux density, averaged 
over the observed frequency range. 
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Fig. 4. a - c; Radial models of chemical abundance as derived by Stauber et al. (2005). The solid curve shows the model of the 
envelope without far-UV and X-ray irradiation by the central object (no-X model), and the dotted curve shows the X-ray model 
including an X-ray luminosity of 8 ■ 10^' erg s~^ and protosteUar far-UV irradiation of Go = 10 at 200 AU fitting best to the 
single-dish data on 24 molecules. Fig. 4d: Density and temperature model of AFGL 2591 used for chemical modeUng (see text). 



sumed to be mostly frozen out on grains, and the total gas-phase 
abundance of sulfur is only 3 • 10"^, maintained by cosmic ray 
sputtering. The abundances of these molecules will thus show a 
jump at the sulfur evaporation radius from a high value inside to 
a very low level farther out. Chemical modeling predicts that the 
abundances of CS and SO in the gas phase are enhanced by X- 
rays after sulfur evaporation. According to the model of Doty et 
al. (2002), the r = 100 K surface has a radius of 1300 AU. Doty 
et al. (2006), on the other hand, predict the water evaporation 
edge at about 670 AU, assuming their preferred age of 8 • 10^ yr 
since YSO formation and using a dynamical model. The mag- 
nitude of the jump depends strongly on X-ray irradiation. In the 
X-ray model the jump amounts to two orders of magnitude for 
CS and to three orders of magnitude for SO. 

HCN is little affected by sulfur and water evaporation. Its 
abundance is greatly enhanced, however, in the central region 
(Sect. 1 and Fig. 4c) due to the presence of atomic N and C, 
and the lack of atomic O. This situation arises through the fol- 
lowing processes: At temperatures above 230 K, the reaction 
with H2 drives most atomic oxygen into OH and finally H2O 
(Charnley 1997; Rodgers & Charnley 2001). On the other hand, 
atomic nitrogen is available in the central region from the dis- 
sociation of N2 by far-UV photons. The association of N by the 
reaction NO + N ^ N2 -1- O is suppressed as all oxygen resides 
in H2O. Thus atomic N and C are available and form CN. As 
CN + H2 ^ HCN + H has an activation energy of 820 K, CN is 
hydrogenated to HCN only at high temperatures. In regions not 
irradiated by X-rays and far-UV photons, atomic carbon forms 
through photodissociation of CO induced by cosmic rays. If pro- 
tosteUar far-UV or X-ray irradiation enhance the C abundance, 
HCN is formed more numerously. Therefore, the HCN abun- 
dance is greatly enhanced by irradiation at high temperatures as 



modeled in the central region within a radius of 250 AU. In out- 
flows or outflow walls direct far-UV irradiation and heating may 
enhance HCN. This is a possible cause of the secondary peaks. 

To find the visibility amplitudes of the models as observed 
by the SMA, the molecular density is first determined by folding 
the fractional abundance with the density model of H2 (Fig. 4d). 
The ID Monte Carlo code of Hogerheijde & van der Tak (2000) 
and molecular data from the Leiden atomic and molecular data 
base (Schoier et al. 2005) are then used to determine the molec- 
ular level populations. Then, radiative transfer is calculated. The 
result is integrated along each line of sight into a two dimen- 
sional map (Fig. 5). Optical depth efl'ects are included, but not 
non-radial inhomogeneity or velocity flows. The map is finally 
transformed into visibility amplitudes applying the (u, v) cover- 
age of the actual SMA observations with the help of MIRIAD. 

The modeling of ro- vibrational HCN(V2=1) emission re- 
quires some assumptions. The level populations were approx- 
imated from the rotational 4-3 transition and weighted by the 
density. 



n(HCN(v2 = 1)) : 



X exp 



n(HCN)g„(HCN(y2 = 1)) 

5«(HCN) 
-[£„(HCN(y2 = 1)) 



kT 



£„(HCN)] j 



(3) 



assuming that the levels are thermalized and in LTE. Thus the 
levels V2 = and 1 have similar abundances at high tempera- 
tures, and the ratio is close to zero at T < 20 K. 

To discuss the effects of abundance variations in radius on 
visibility amplitudes, we recall four basic properties of Fourier 
transformation. (/) Adding a small quasi-constant component 
to the abundance enhances just the zero value and leaves un- 
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Fig. 5. Line maps of chemical model including X-rays and far- 
UV irradiation as if observed by the SMA. Contours are at 10% 
steps of the peak. The 50% curve is bold. The SMA beam is 
shown for comparison in Fig. 5a. 



changed the visibilities at non-zero baselines, (ii) The visibility 
amplitudes of a partially resolved peak depend linearly on the 
flux of this source. (Hi) Expanding a resolved Gaussian source in 
size reduces the non-zero amplitudes. A jump in the density acts 
like an additional source having the size of the jump radius and 
enhances the visibilities below the corresponding baseline (Eq. 
2). (/v) A boxcar density model transforms into a s/nc-function 
in visibility ampUtude. A boxcar molecular abundance distribu- 
tion out to a radius r thus causes nulls in visibility amplitude at 
xu given by Eq. (2). These properties of Fourier transformation 
persist basically in the case of a power-law given by the den- 
sity model. An abundance jump at 1300 AU predicted by the 
model for CS and SO thus yields a zero in visibility ampUtude 
at 113 kZ. The value can differ slightly if the density model does 
not have an exact boxcar shape. Visibilities observed for AFGL 
2591 in the range of 38 to 215 kA are sensitive indicators of spa- 
tial structures with radii from a few hundred to a few thousand 
AU. 



6. Comparison of models with data 

There are four ways to compare the chemical model of Stauber 
et al. (2005) with the data: 

- The total flux of the model in a 14" x 14" area can be read- 
ily compared with the single-dish observations by the JCMT 
(both given in Tab. 2). The total JCMT flux corresponds to 
the value at zero kA in Fig. 3 and includes in particular the 
radii 3500 - 7000 AU, missing in the SMA total flux. As the 
chemical model was fitted to 24 molecular column depths 
derived from single-dish data including the lines selected 
here, it is not surprising to find here general agreement. The 
largest deviation is the main isotope HCN, observed by the 
JCMT more than five times brighter than predicted. 



- The region with radii 600 - 3500 AU is best compared in 
visibiUty ampUtudes (Fig. 3). It is of special interest as it 
contains the region influenced by protostellar X-ray emission 
as well as the sulfur evaporation region. 

- The center region (< 600 AU) is related to the largest (m, v) 
distances (200 - 220 kA, Tab. 2). Alternatively, the 0.5" cen- 
tral field of the model can be tested by the central SMA beam 
measurement (marked 0.5" in Tab. 2). The flux observed in 
the SMA beam at the peak of the continuum agrees with the 
unresolved flux within the given accuracy except for the flat 
SO emission (see next paragraph). 

- Finally, the predicted size of the model depicted in Fig. 5 can 
be contrasted to the observed map (Fig. 1) and the Gauss fit 
(Tab. 1). 

A first, important result: The molecules reacting to sulfur 
evaporation, CS and SO, have larger emitting regions than the 
HCN lines at the 50% isophote in Fig. 1 as expected from the 
X-ray model (Fig. 5). CS has a smaller diameter measured by 
the Gauss fit than SO (Tab. 1). In effect, the CS emission con- 
tains less flux in the central peak and has a higher peak. Thus its 
shape is spikier than the SO emission. CS appears clearly more 
extended than HCN. In the foUowing we distinguish between the 
size at the FWHM level and that of a broad pedestal, and discuss 
the lines individually. 

6.1. CS line emission 

In general, the modeled CS visibilities agree with the observed 
ones. There are three deviations: (i) The kA value is a factor 
of 3 below the observed one (JCMT flux. Tab. 2). This indicates 
that the model underestimates the flux at radii 3500 - 7000 AU. 

The CS model produces marked minima at 105 and 190 kA 
where the amplitude reaches nearly zero (Fig. 3a). It is the re- 
sult of the extremely sharp abundance jump by two orders of 
magnitude at the sulfur evaporation edge (X-ray model. Fig. 4a), 
causing a null in Fourier space at 102 kA and its approximate 
integer multiples. ( Hi) Figure 5 predicts a FWHM size of about 
2.1", remarkably larger than 1.04", as measured from the main 
emission peak (Tab. 1) convolved with the beam. 

6.2. SO line emission 

The model for SO 88-7? emission also agrees relatively well with 
the data. In particular, we note that the kA value agrees with the 
observed one (JCMT flux. Tab. 2). Most importantly SO being 
our best X-ray tracer, the X-ray model reproduces the observed 
data much better than the no-X model. The SO abundance de- 
creases below 300 AU (Fig. 4b), reducing the optical thickness in 
the central region and diminishing the Fourier power of the un- 
resolved central region. Thus the SO visibilities at >200 kA are 
predicted to be relatively low and less peaked than CS. Indeed, 
the observations show no signature of an unresolved central peak 
contrary to CS. 

There are two major deviations: (;) Like for CS, the model 
impUes a jump in SO abundance at 1300 AU by more than a 
factor of 1000 (Fig. 4b). The deep minima in Fourier space are 
not observed. Nevertheless, a dip in the observed visibilities at 
1 15 kA is suggestively close to the predicted 1 17 kA (Fig. 3b). If 
real, it may be interpreted by a jump in abundance. The dip may 
be visible only in SO where the jump is predicted to be larger 
than in CS. (ii) On the other hand, the predicted FWHM size is 
1.9" (Fig. 5), a factor of about 1.6 larger than observed (Tab. 1, 
Fig. 1). If the size were determined by a jump, the dip would be 
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expected at 1.6 times smaller (m, v) distance, again contrary to 
observations (Fig. 3). 

6.3. HCN line emissions 

The main isotope HCN 4-3 model underestimates the total flux 
at the 7000 AU radius scale by a factor of 5. Thus HCN has an 
additional intense pedestal emission collected by the JCMT sin- 
gle dish, but not predicted by the model. On the other hand, the 
model fits roughly in the inner part (< 1 100 AU, corresponding 
to > 120 k/l). The predicted FWHM of 0.8" (Fig. 5) agrees with 
observations. 

Contrary to the rotational HCN fine, the modeled ro- 
vibrational HCN(v2=l) total flux agrees well with the JCMT ob- 
servations. However, the visibility amplitudes over the range 60 
to beyond 210 k/l (radii 600 - 2200 AU) are overestimated, sug- 
gesting a main peak size larger than modeled by about a factor 
of 3. This is corroborated by the flux in the central SMA beam, 
being smaller than predicted (Tab. 2). 

The HC'^N model follows the main isotope HCN in many 
properties. It fits relatively well with the data in the inner part, 
predicting even the unresolved peak fiux to within better than 
a factor of two. The model fits well at large visibiUties (Fig. 
3, » 600 AU), and in the unresolved peak. There is also good 
agreement with the observed FWHM size. However, the same 
shortcoming persists: The total flux is a factor of 3 lower than 
observed by the JCMT single-dish. It requires an additional ex- 
tended component at radii > 1 100 AU (< 120 k A). 

In conclusion, HCN and HC'^N agree qualitatively with 
the predicted central abundance peak, but require an additional 
broad pedestal emission beyond about 1100 AU. This broad 
component is not clearly perceptible in HCN(v2=l) visibilities, 
although the secondary emission peak is detected also in this line 
(Fig. 1). 



7. Discussion 

7.1. X-ray and far-UV irradiation 

The question of X-ray emission by the deeply embedded HMPO 
AFGL 2591 can be addressed by these SMA observations in 
more depth than previous studies based only on single-dish ob- 
servations. Spatial resolution allows testing independently if and 
how X-rays affect molecular abundances with distance from the 
protostar, depending on density, temperature and irradiation. The 
X-ray model of Stiluber et al. (2005) includes also far-UV irra- 
diation, which has Uttle effect beyond 300 AU in a spherically 
symmetric model. If present. X-ray and far-UV irradiation have 
similar chemical effects, differing in the primary ionization and 
in penetration depth. In the inner part (< 1000 AU) of a spher- 
icaUy symmetric model, SO mainly traces X-rays, and HCN 
mainly traces far-UV irradiation (Sect. 1). CS can be enhanced 
by both irradiations (Stauber et al. 2004, 2005). 

These SMA observations support the hypothesis of X-ray 
irradiation in two ways: Firstly, Fig. 3b demonstrates that the 
model including X-rays fits the SO visibility amplitude better 
than the no-X model. The difference is remarkable at large radii 
( > 1300 AU, corresponding to < 100 kA) where X-rays com- 
pletely dominate far-UV irradiation. Secondly, the CS emission 
is predicted to have a large main source in the presence of X- 
rays. Figure 4a shows that this is the result of X-ray irradiation 
within the sulfur and water evaporation region and does not oc- 
cur without X-rays. At radii larger than about 1000 AU, the CS 



source size becomes therefore a further X-ray tracer. The mea- 
surements are consistent with this prediction and thus corrobo- 
rate X-rays. 

The HCN abundance is sensitive to far-UV irradiation, more 
than to X-rays (Section 1). The HCN lines deviate from the mod- 
eled visibility amplitudes mostly in the 500 - 3500 AU radii 
range (Fig. 3) at distances comparable to the secondary emis- 
sion peak possibly located in the outflow. In addition, HCN and 
HC'^N, but not HCN(v2=l), contain a broad component beyond 
3500 AU accounting for nearly 90% of the main isotope HCN 
emission. The HCN(v2=l) line is preferentially emitted in the 
very irmer part. This appears to be the result of its excitation 
mechanism. It is supported by the high brightness temperature of 
the HCN(v2=l) line (Tab. 1), consistent with observed HCN ab- 
sorption lines in the infrared indicating temperatures up to 1300 
K (Boonman et al. 2001) and its small source size. 

Far-UV radiation propagates only in a low density gas, but 
scatters more easily than X-rays. The deviations of HCN and 
the CS from the predictions of a ID envelope model may be 
explained as additional emission in outflows or outflow walls. 
HCN emission in high-mass outflows has recently been reported 
by Zhang et al. (2007). In addition to far-UV irradiation, HCN 
at large radii requires higher densities and temperatures than the 
model values. The secondary emission peak (Fig. 1) contributes 
to this pedestal emission. The far-UV origin of the emission 
is further supported by the absence of SO emission at the sec- 
ondary peak, as SO is not much enhanced by far-UV. 

Shocks may be an alternative cause for secondary peaks and 
pedestal emission. The hypothesis is supported by the observed 
line asymmetry in main isotope HCN. However, the other lines 
being relatively narrow indicate that most of the structure ob- 
served by the interferometer does not originate from shocks. 
Moreover, single-dish observations do not indicate shocked gas 
(e.g. van der Tak et al. 2003). 

7.2. Sulfur evaporation 

In addition to protosteUar irradiation, the SMA observations also 
probe sulfur evaporation. The observations indicate that CS and 
SO line emission extend well into the 900 - 2500 AU radii re- 
gion (Figs, la and b, Tab. 1). A pronounced abundance enhance- 
ment in the range < 3500 AU reduces the ratio of central beam to 
total SMA flux. Thus, the differences between CS and SO com- 
pared to HCN in Tab. 2 are consistent with the CS and SO cen- 
tral peak sizes being more extended due to sulfur evaporation. 
Such a feature is less pronounced or absent in all three HCN 
fines. The amount of the abundance jump of CS and SO in the 
chemical model of Stauber et al. (2005), however, does not agree 
when compared with observations in terms of visibilities. The 
absence of pronounced minima in CS and SO visibilities (Fig. 
3) may require a smaller jump than predicted. Instead of reduc- 
ing the jump magnitude, the effects of the abundance jump may 
be smoothed by inhomogeneity. If the sulfur evaporation edge 
varied in radius by +10% for different longitudes and latitudes, 
the superposition of the difl'erent visibility distributions would 
suffice to eliminate the spurious minimum. 

The predicted source sizes of CS and SO are largely deter- 
mined by the position of the evaporation edge in radius. They 
exceed the observed sizes of the main peak considerably. Could 
this be caused by a sulfur evaporation temperature different from 
the assumed value of 100 K? To quaUtatively explore the effect, 
the chemical model was modified to allow sulfur evaporation at 
60 K instead of 100 K. Such a change would double the sulfur 
evaporation radius and augment the region of enhanced CS and 
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SO abundances. The more distant evaporation edge, however, 
would make the source even larger than observed and would 
move the nulls in the model visibilities in Fig. 3a and 3b to (u, v) 
distances smaller by a factor of two. This would diminish the 
quality of the fit. The SMA observations suggest a sulfur evap- 
oration radius in the range 700 - 1300 AU. A higher evapora- 
tion temperature is unlikely in view of laboratory experiments 
(Section 1). Thus, the observations are in agreement with sulfur 
evaporation near 100 K, roughly coinciding with water evapora- 
tion. 

The CS and SO main peak may appear smaller in size be- 
cause of the above spatial smoothing or a different temperature 
model. Note that a distance of AFGL 2591 larger than the as- 
sumed 1 kpc does not solve the discrepancy, as the angular size 
of the 100 K radius is independent of distance (S.D. Doty, per- 
sonal communication). 

7.3. The innermost envelope 

The agreement of the model with SMA observations of the main 
peak is acceptable for the unresolved part of the star forming 
region (< 600 AU in radius). It emits much less than the re- 
solved region, for example less than 5% of the total main iso- 
tope HCN emission. The chemical models overpredict the unre- 
solved HCN(v2=l) by a factor of three, and underestimate the 
unresolved part of the main isotope HCN. The observed strong 
central peak of the HCN conforms with infrared observations by 
Boonman et al. (2001), indicating an abundance of 10"^ (Sect. 
1). 

We have estimated the optical thickness of the lines from the 
model density (taken from van der Tak et al. 1999) and abun- 
dance (from Stauber et al. 2005). The column density becomes 
insensitive to temperature above about 100 K for lines with low 
upper state energy (thus excluding HCN(v2=l); van der Tak et 
al. 2006, Fig. 5). The radiative transfer in molecular lines was 
calculated applying non-LTE excitation and the escape probabil- 
ity method (RADE^fl van der Tak et al. 2007). All the observed 
lines, even HC'^N, become optically thick in the center, as noted 
independently in Fig. 2 for some of the lines. Because of opti- 
cal thickness and since some of the emission is resolved out by 
the interferometer, we refrain from estimating molecular column 
densities based on the SMA observations. 

8. Conclusions 

Observations of AFGL 2591 at subarcsecond resolution have 
been made using the SMA in the frequency range 343 - 354 
GHz. The resulting baselines are ideal to explore the envelope of 
the high-mass YSO at 600 - 3500 AU radius, a region irradiated 
by protostellar X-rays and far-UV photons. The observed radii 
extending over a temperature range of 50 - 150 K include also 
the radius of sulfur evaporation from grains. Chemical models 
predict large changes in abundance of some molecules, such as 
CS and SO, in this range. The magnitude of the change strongly 
depends on the flux of possible X-ray iiTadiation (Fig. 4). The 
main results are: 

- The visibility amplitudes of SO, the best X-ray tracer ac- 
cording to the model of Stauber et al. (2005), matches the 
observations much better when X-rays are included (Fig. 
3b). With the exception of HCN(v2=l) at small radii (* 600 
AU), this holds as well for all other observed molecules. 

^ http://www.sron.rug.nl/~vdtak/radex/radex.php 



According to the chemical models, the CS abundance is en- 
hanced in the inner envelope and has a significant jump at the 
sulfur evaporation edge only under the influence of consid- 
erable X-ray irradiation (Fig. 4a). Thus the detection of sig- 
nificant extended CS emission in the main peak also corrob- 
orates X-ray irradiation. Therefore, these SMA observations 
support the existence of protostellar X-rays in this HMPO. 

- However, the X-ray model of Stauber et al. (2005) predicts 
pronounced minima in visibility amplitude at a (u, v) dis- 
tance near 110 kA for CS and SO that are not observed. 
Nevertheless, SO, predicted to have the larger jump, shows a 
dip at about this baseline that tentatively may be interpreted 
as a weak signature of it. The model still has to be modi- 
fied considerably to agree with the CS and SO observations. 
The jump may be reduced by inhomogeneity in temperature 
at the radius of sulfur evaporation, varying the radius in a 
range of at least +10%, so that the minimum is averaged out 
and the apparent source size increases. Evolutionary effects 
are alternatives, producing a time delay of grain heating and 
sulfur evaporation in a non-stationary model. As water evap- 
orates at a similar temperature and radius, a reduction of the 
jump may then also be expected for H2O. 

- The size of the CS and SO main sources indicate an increase 
of CS and SO abundance in the range of 100 - 190 kA, cor- 
responding to 1300 - 700 AU. The results thus concur with 
sulfur evaporation at 100 K, similar to desorption of water 

- The relatively small ratios of central beam SMA flux to to- 
tal SMA flux (Tab. 2) confirm that CS and SO indeed extend 
beyond 500 AU. Since the observed FWHM size of the emis- 
sion is smaller than predicted, this suggests a sulfur evapora- 
tion edge at a smaller radius ( < 1000 AU) than modeled. 

- On the other hand, the visibility amplitudes of CS, HCN, and 
HC'^N (but not SO and HCN(v2=l)) drop in observed vis- 
ibility from kA (single-dish value) to the lowest (u, v) dis- 
tance (30 - 40 kA) more than predicted (Tab. 2), requiring in 
addition a broader source component, or pedestal, reaching 
up to 7000 AU at a low level. 

- A secondary emission peak located 1 180+140 AU northwest 
of the main peak within the approaching outflow cone is part 
of a pedestal emission of the HCN lines from the outflow or 
its wall. The secondary peak is seen in all HCN lines and 
also in CS, a molecule that has been found to be sensitive 
also to far-UV in chemical modeling. The secondary peak 
is not observed in SO which is insensitive to far-UV irradi- 
ation. The pedestal emission is clearly indicated also in the 
visibility amplitudes of CS, HCN, and HC'^N, but not SO 
and HCN(V2=1). 

- The line emission from the outflow or its walls is in all cases 
less than 30% of the total interferometric flux. 

- The visibility analysis of the SMA data indicates that CS 
emission has a more pronounced unresolved central peak 
than SO. This agrees with the predictions of the chemical 
models (Figs. 4a and b). The observed unresolved emission 
is slightly larger than modeled for CS and the main isotope 
HCN. 

In conclusion, we present the first quantitative spatially re- 
solved tests of an X-ray model for YSO envelopes previously 
fitted to unresolved single-dish data (Stauber et al. 2005). The 
SMA observations support the hypothesis of protostellar X-rays 
and far-UV irradiation, but require an additional broad pedestal 
component possibly located in the outflows. Overall, these ob- 
servations demonstrate the power of submillimeter interferome- 
ters for directly testing abundance profiles of molecules in pro- 
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tostellar envelopes on scales where many important physical and 
chemical processes are happening. 

For a quantitative description of the radial abundance varia- 
tions in the envelope it will be necessary to fit a detailed, possibly 
non-spherically symmetric model to imaged data. This would re- 
quire a large computational effort, beyond the scope of this work, 
but could potentially provide a more detailed picture for the in- 
fluence of X-ray and far-UV protostellar emission compared to 
other chemical effects in high-mass protostars. 
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